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The ability to distinguish between high and low levels of task engagement in the 
real world is important for detecting and preventing performance decrements during 
safety-critical operational tasks. We therefore investigated whether functional Near 
Infrared Spectroscopy (fNIRS), a portable brain neuroimaging technique, can be used to 
distinguish between high and low levels of task engagement during the performance of 
a selective attention task. A group of participants performed the multi-source interference 
task (MSIT) while we recorded brain activity with fNIRS from two brain regions. One 
was a key region of the "task-positive" network, which is associated with relatively high 
levels of task engagement. The second was a key region of the "task-negative" network, 
which is associated with relatively low levels of task engagement (e.g., resting and not 
performing a task). Using activity in these regions as inputs to a multivariate pattern 
classifier, we were able to predict above chance levels whether participants were engaged 
in performing the MSIT or resting. We were also able to replicate prior findings from 
functional magnetic resonance imaging (fMRI) indicating that activity in task-positive and 
task-negative regions is negatively correlated during task performance. Finally, data from a 
companion fMRI study verified our assumptions about the sources of brain activity in the 
fNIRS experiment and established an upper bound on classification accuracy in our task. 
Together, our findings suggest that fNIRS could prove quite useful for monitoring cognitive 
state in real-world settings. 
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INTRODUCTION 

The ability to distinguish between high and low levels of task 
engagement is important for detecting and preventing perfor- 
mance decrements during safety-critical operational tasks in the 
real world. Examples of such tasks include commercial aviation, 
monitoring for air traffic control, executing space walks, perform- 
ing surgery, and driving. Since accident-causing errors can be 
made even by skilled professionals (Dismukes et al., 2007), the 
ability to monitor cognitive state measures for low levels of task 
engagement in real time could be useful for developing an "early 
warning system" for detecting and preventing performance errors 
before they occur. 

The use of cognitive state measures to optimize human per- 
formance (for example by informing flight automation or the 
operator themselves of a hazardous state) has been of particu- 
lar importance to aviation safety (Pope et al., 1995; Schnell et al., 
2004) and to the Augmented Cognition program (Raley et al., 
2004) of the Defense Advanced Research Projects Agency of the 
United States. Additionally, such research is highly relevant to 
space flight, since adaptation to microgravity can cause perfor- 
mance decrements due to motion sickness, lack of sleep, loss of 
sensorimotor control, increased stress or mood changes (Cowings 
et al, 2003). More generally, the ability to monitor cognitive state 
for low levels of task engagement could be helpful for detect- 
ing and preventing vigilance decrements due to sleep-deprivation 



(Drummond et al, 2005; De Havas et al., 2012) or distraction 
(Strayer etal.,2011). 

Monitoring brain activity may provide an effective means for 
monitoring cognitive state, and in particular for distinguishing 
between high and low levels of task engagement. Numerous func- 
tional magnetic resonance imaging (fMRI) studies have revealed 
that activity increases in a so-called "task-positive" network, 
which includes the dorsolateral prefrontal cortex (DLPFC), dorsal 
anterior cingulate cortex (dACC), superior and inferior parietal 
lobe (SPL and IPL), and anterior insula (Al), when participants 
perform a task as compared to when they rest (MacDonald et al., 
2000; McKiernan et al, 2003; Dosenbach et al, 2006). In con- 
trast, activity increases in a so-called "task-negative" network, 
which includes the anterior medial frontal gyrus (aMFG), poste- 
rior cingulate cortex (PCC), and certain regions of lateral parietal 
cortex (LPC), when participants rest as compared to perform a 
task (Raichle et al., 2001; Greicius et al, 2003). In other words, 
activity in the "task-positive" and "task-negative" networks is 
negatively correlated (Fox et al., 2005; Kelly et al., 2008). For this 
reason, monitoring activity in key regions of the "task-positive" 
network alone or monitoring activity in key regions of the "task- 
positive" and "task-negative" networks together may distinguish 
between relatively high and relatively low levels of task engage- 
ment (Drummond et al, 2005; Weissman et al., 2006; Chee et al., 
2008). Given recent data indicating that interactions between key 
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regions of the "task-positive" and "task-negative" networks vary 
with task engagement (Prado and Weissman, 201 1), we predicted 
that either approach for monitoring brain activity would allow us 
to distinguish between relatively high and relatively low levels of 
task engagement, but that the latter approach would likely prove 
most effective. 

Since hemodynamic activity cannot be monitored with fMRI 
outside of a laboratory, we employed functional Near Infrared 
Spectroscopy (fNIRS) to determine whether monitoring brain 
activity is an effective method for monitoring cognitive state. 
FNIRS is a portable optical neuroimaging technique that can 
be used to quantify hemodynamic activations. Moreover, it is 
relatively low-cost, non-confining, non-invasive, and safe for 
long-term monitoring (Boas et al, 2004; Gibson et al., 2005; 
Gratton et al, 2005; Steinbrink et al, 2005; Schroeter et al, 2006). 
Finally, temporal resolution is sub-second, and spatial resolution 
is on the order of 1 cm 2 at best (Strangman et al., 2002a; Obrig 
and Villringer, 2003; Bunce et al., 2006). Measurements have been 
shown to be consistent with fMRI (Kleinschmidt et al, 1996; 
Strangman et al., 2002b; Steinbrink et al., 2005; Huppert et al., 
2006; Schroeter et al, 2006; Emir et al., 2008) and electroen- 
cephalography (EEG) (Moosmann et al, 2003; Li et al, 2007). 
Critically, the ambulatory nature of fNIRS allows neu- 
roimaging in the field. Thus, fNIRS may evolve into a 
synergistic complement to EEG and other physiological 
measures for monitoring cognitive state during operational 
tasks. 

In the present study, we employed fNIRS to determine whether 
it is possible to distinguish between high and low levels of task 
engagement (i.e., performing a task vs. resting). Specifically, we 
monitored brain activity from the DLPFC in the "task-positive" 
network and from the MFG in the "task-negative" network while 
participants alternated between performing and not performing 
a cognitive task. We then employed multivariate pattern classifi- 
cation techniques in an effort to distinguish between periods of 
task performance and periods of rest. 

To facilitate our ability to make this distinction, we asked par- 
ticipants to perform the multi-source interference task (MSIT; 
Stins et al, 2005; Bush and Shin, 2006). The MSIT is a selec- 
tive attention task in which optimal performance requires par- 
ticipants to suppress multiple sources of interference (Stroop, 
Eriksen, and Simon). Thus, it reliably and robustly activates the 
"task-positive" network, even in individual task blocks from the 
same participant (Bush and Shin, 2006). Given these character- 
istics, we reasoned that the MSIT would provide a strong signal 
with which to monitor task engagement. Consistent with this 
reasoning, in the present study we were successful at distinguish- 
ing between relatively high and relatively low periods of task 
engagement. 

We also conducted a companion fMRI study for verification 
and comparison purposes. First, given that DLPFC and MFG 
activity was recorded with fNIRS at the scalp surface, we wished 
to verify that our paradigm actually elicited hemodynamic acti- 
vations in these regions. Second, we wished to compare fNIRS 
with fMRI with regard to the ability to distinguish between high 
and low levels of task engagement. Given that fMRI detects motor 
cortex activation reliably enough for routine use (Moller et al., 



2005), and given that the motor cortex should be activated by 
the button presses required by our task, we expected that includ- 
ing such activation as an input to our classification algorithms 
would produce the highest levels of accuracy. Thus, we reasoned 
that classification accuracy with fMRI when including motor cor- 
tex activation would establish an "upper bound" for classification 
accuracy expectations. As expected, the fMRI study confirmed 
both of these predictions. 

METHODS 
BEHAVIORAL METHODS 

Participants 

Seven participants (three females, four males) completed the 
fMRI study. Five participants (two female, three male) completed 
the fNIRS study (including two from the fMRI study). All partic- 
ipants practiced the behavioral task with performance feedback 
for 1 min at the beginning of the study. Next, they performed 
a set of four 7-min-long runs. Four runs each were completed 
for the fNIRS and fMRI experiments, which were performed on 
separate days. Human participant data were collected according 
to a protocol approved by both the University of Michigan IRB 
MED and the NASA IRB. Informed consent was obtained from 
all participants, who were healthy adults between the ages of 21 
and 50 years. All participants were right-handed but one, who was 
ambidextrous. 

Behavioral task 

In each trial of the MSIT (duration, 2 s), participants viewed 
a horizontally-oriented array of four digits at the center of the 
screen (duration, 1 s). A target digit was printed in a larger font 
than each of three distracter digits (72 point vs. 60 point). The 
stimulus sizes were chosen to make the digits visible within the 
MRI scanner. Moreover, the large target digit was chosen to be 
about 20% larger than the small distracter digits to make the task 
sufficiently difficult. Indeed, Stins et al. (2005) observed a much 
smaller interference effect in the MSIT when the large target digit 
was 33% larger than the small distracter digits. Since the goal 
of the present MRI experiment was to contrast activity during 
task performance to activity during rest, we wanted to ensure that 
the task was difficult enough to elicit a good deal of task-related 
activity. 

Participants were instructed to identify the target digit (1, 2, 3, 
or 4) by pressing a key with one of four fingers ( 1 = right thumb, 
2 = right index finger, 3 = right middle finger, 4 = right ring 
finger) as quickly as possible without making mistakes. For con- 
gruent trials (50%), both the spatial position of the target and 
the identity of the three distracter digits were mapped to the cor- 
rect response (1 111, 2222, 3333, 4444). For incongruent trials 
(50%), both the spatial position of the target and the identity of 
the three distracter digits were mapped to a conflicting response 
(2 1 1 1 , 2 1 22, 3343, 4443 ) . Participants alternated between trials 
composed of Is and 2s and trials composed of 3s and 4s across tri- 
als to prevent immediate stimulus and response repetitions (Mayr 
et al, 2003; Jimenez and Mendez, 2013). 

Responses were recorded via the keyboard (the "n," "u," "9," 
and "0" keys) of a laptop used in fNIRS experiments and via a 
MR-compatible response device in fMRI experiments. The task 
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was implemented using a combination of MATLAB, 2012 (Natick, 
MA) and the Psychophysics Toolbox (Kleiner et al, 2007). 

Functional neuroimaging experimental design 

We employed a block design. In each of four runs, an initial 16 s 
rest block was followed by 1 2 alternations between the MSIT ( 1 6 s 
— 8 trials per block) and rest (16 s). Each run lasted 400 s. 

fNIRS METHODS 
fNIRS data acquisition 

Hemoglobin concentration changes were measured using an 
Imagent NIRS instrument and fiber optic cables (ISS, Inc.). 
Eleven rigidly-connected source-detector pairs were used, and 
each source fiber delivered both 690 and 830 nm wavelength light. 
The data collection rate was 6.25 Hz. During each task run, 2500 
time points were collected. Eight sources were located around 
one detector placed over the DLPFC region, and three sources 
were located around a second detector placed over the MFG. The 
array of head probes and fiber optic cables used to interrogate 
the MFG in this study is shown in Figure 1, right. The sources 
were held in place using both clear and blacked-out plastic at 
the locations shown with respect to the International 10-20 loca- 
tions in Figure 1, left. The array of probes used to interrogate 
the MFG contained two sources placed 3 cm from the detector 
placed between FPz and FP2. The array of probes used to inter- 
rogate the DLPFC contained seven sources placed 3 cm from the 
detector placed near F4. The clear material (not shown) improved 
visual inspection for placement and hair control at the probe- 
skin interface. Blackout material was applied over and around the 
probes to block ambient light. Motor regions were not simultane- 
ously interrogated due to instrumentation and fiber optic probe 
limitations. 

The probes were located with the aid of an electroencephalog- 
raphy net (64-channel HydroCel Geodesic Sensor Net by EGI, 



Inc.) applied according to EGI Inc.'s instruction. The nets were 
used to identify the International 10-20 locations for each par- 
ticipant at each visit, not to record EEG data. A mark was made 
under the same net pedestals for all participants. The use of the 
nets made localization reliable, consistent and expeditious. The 
MFG mark was placed under a pedestal half way between FPz 
and FP2. This placed the sources for the MFG array (which were 
2 cm from each other) at about the midline, with the detector 
then about 3 cm from the midline, to avoid the superior sagit- 
tal sinus as much as possible. The DLPFC mark was placed at the 
pedestal for EEG channel 59, which is immediately inferior to F4. 
The DLPFC array was placed by hand such that the detector and 
the source for channel 3 straddled its mark. The MFG array was 
placed similarly for channel 2. The probe arrays were secured with 
Velcro straps. The probes were not moved between the four runs 
unless the participant requested adjustment for the purposes of 
comfort. In those cases, care was taken to relieve pressure on the 
head without translating the probes. 

Each probe array included a shallow source located 1 cm from 
the detector for superficial physiological and nuisance signal 
regression. The sources located 3 cm from the detectors pro- 
vided deep traces of interest which sampled brain tissue, while 
the sources at 1 cm provided shallow traces which sampled only 
superficial tissue due to the proximity to the detector (Gagnon 
et al, 2011). The shallow channels were primarily sensitive to 
physiological changes in the skin, while also being sensitive to 
nuisance signal contributions such as motion of the rigidly- 
connected probe and ambient light exposure. 

The Imagent instrument employs photomultiplier tubes for 
optical intensity detection. Gain settings for these were set on a 
per-participant basis after applying the probes and ensuring most 
of the hair was parted under the probe tips. Gain was increased 
and probes were re-adjusted to make better contact with the skin, 
iteratively, until as many channels as possible detected continuous 
wave intensity signal above a threshold of 500 analog-to-digital 
counts. None of the channels with low signal produced the best 
fit to the task model, and thus were not passed to the classifi- 
cation step (described below). The shallow channel sources at 
1 cm may damage the detectors at gain settings appropriate for 
sources at 3 cm due to the lack of signal attenuation over their rel- 
atively shorter optical path length. To avoid this, delivered optical 
intensity was reduced for the shallow channel sources by layer- 
ing optically absorbent pigment on partially-transmissive tape 
between the ends of the fiber probe tips and the skin. 

fNIRS data processing 

Both oxygenated and deoxygenated hemoglobin concentration 
([Hb]) changes were calculated from the filtered raw contin- 
uous wave intensity measured traces using the Modified Beer 
Lambert Law (Delpy et al., 1998), then normalized (Huppert 
et al, 2009). Filtering was set to include 0.008-0.08 Hz to focus 
on sustained task activations while removing very slow drift and 
higher frequency physiological and motion contributions. 

Standard linear regression was then used to remove physio- 
logical contributions from the measured signal to produce the 
functional task signal, and to select probes from each array. Probe 
channels, and specifically the oxygenated or deoxygenated [Hb] 
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FIGURE 1 I Head probes used to interrogate the DLPFC and MFG. Left: 

source channels (red) are arranged at a separation distance of 3 cm from 
two detectors (green, "D"), shown schematically with respect to the 
International 10-20 locations. The shallow sources, located at 1 cm from 
the detectors, provide channel 4 for the DLPFC array, which is near F4, and 
channel 2 for the MFG array, which is between FPz and FP2. The probes 
consistently producing the traces with the highest task model fit 
parameters were DLPFC probe channels 1-3 and MFG channel 1 (as 
numbered, left). Right: the skin-side of the array of head probes used for 
the MFG. 
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FIGURE 2 | Classification accuracy averaged across four runs for each 
of five participants. Two fNIRS time traces were used as support vector 
machine input features in each case. Left-hand bars: within-network pairs 
(green). Right-hand bars: cross-region pairs (red). Accuracy is the number 
of time points for which the prediction matched the truth label out of all 
2500 time points. An accuracy of 50% represents prediction by chance. 
Error bars represent ± one standard deviation. 
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FIGURE 3 | Two across-region time traces averaged across all 20 runs. 

These traces illustrate filtered, normalized, and corrected (see fNIRS 
methods) hemoglobin concentration changes. For the purposes of this 
group average, all deoxygenated traces were first inverted. Truth labels are 
indicated by the green trace with task at +1 and rest at —1. The 
across-network correlation for the group was —0.58. Task blocks were 16s 
each, with equal rest time. Dotted lines show ± one standard deviation for 
each time point. 



trace from that probe, were selected for use in classification based 
on how well the measured deep traces fit the expected functional 
activation task response, quantified by their beta fit parameter. 
The expected task response was modeled by convolving the boxcar 
task onset signal with the hemodynamic response function. The 
expected deoxygenated [Hb] time series was set to the negative 
of the oxygenated [Hb]. The shallow trace was smoothed using 
a 6-timepoint moving average, and task-like response in it was 
removed (by a separate regression step) prior to use as a nuisance 
regressor in the design matrix. The shallow trace removal was per- 
formed within-species. That is, the oxygenated [Hb] trace from 
the shallow channel was regressed from the measured oxygenated 
trace of each deep channel on that detector's array, and similarly 
for the deoxygenated traces. The functional task signal, which was 
effectively the measured signal minus the fitted nuisance signal, 
was then passed to the classification step. 

Within- and across-network regional correlations, respec- 
tively, were defined as the time-series correlation coefficients 
relating activity between different regions of the DLPFC in the 
task-positive network and between the DLPFC and the MFG in 
the task-negative network. Each of these correlation values was 
averaged across all four runs in each participant, before it was 
averaged across participants. All statistical tests on correlation 
and accuracy values were one-tailed, with comparisons paired by 
participant. 

fNIRS classification 

Classification was performed using two traces as input features 
to Support Vector Machines (SVM). Scripts were implemented in 
MATLAB, 2012 using LibSVM (Chang and Lin, 2011). We choose 
SVM for good performance with ease of implementation, pro- 
cessing speed in the interest of future real-time application, and 
the ability to use tuning parameters to optimize feature separation 
for each participant. First, the best two traces for each participant 
were selected from the array of seven deep DLPFC channels (a 
within-network pair, noted as DLPFC and DLPFC2 in Figure 2), 
based on which traces (oxygenated or deoxygenated) best fit the 
functional activation model as described above (which had the 
highest beta fit parameter considering all four runs). The traces 
which best fit the task model for each participant were assumed 
to make the best input features for producing the highest classi- 
fication accuracy. The same best DLPFC trace was then paired 
with the best trace for each participant from the array of two 
deep MFG probes (an across-network pair, noted as DLPFC and 
MFG in Figures 2, 3) and a separate, second classification step 
was performed. Trace selections were not changed across runs. 

A SVM model was trained to discriminate high from low lev- 
els of task engagement using three of each participant's four runs. 
Its prediction accuracy was tested on the participant's fourth run. 
Thus, training and prediction were always conducted within par- 
ticipants. All permutations were computed for each participant, 
such that the SVM model's prediction accuracy could be deter- 
mined for each of the four runs. As described above, this was 
done separately for within- and across-network input pairs. The 
truth labels used for training and accuracy determination pur- 
poses were determined from the boxcar task onset signal, but 
shifted 4 s later to account for the delay of the hemodynamic 
response (the green trace in Figure 3). 



A radial basis function kernel was employed for the SVM, 
which involves two tuning parameters. Kernel parameter gamma 
(g) determines the non-linearity of the mapping of our few fea- 
tures into a multidimensional space in which the prediction class 
labels are determined. A lower cost of error (c) allows more flex- 
ibility and generalizability of the SVM model by reducing the 
cost of misclassifications during model training (Chang and Lin, 
2011). Non-linear mapping may improve performance depend- 
ing on g. To determine the best possible prediction accuracy 
achievable with the methods of this study, c was tested at 0.1, 
1, and 10, and g was tested at 0.001, 0.01, and 0.1 (a total 
of nine cases) in each classification step to optimize accuracy. 
These values were selected after running a test case varying c 
and g each across nine orders of magnitude. Accuracy did not 
change appreciably across four orders of magnitude, and the c, 
g values used were selected from this region. The highest clas- 
sification accuracies produced using the same set of parameters 
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across all four runs (the same c and g, given in Table 1) were 
selected for each participant and reported (see Figure 3). The 
best c, g parameters found upon optimization are not extreme 
in value and are not identical across participants. Accuracy was 
determined by comparing the work or rest state predicted for 
each time point to the known states (i.e., performing the task 
or resting). 

fMRI METHODS 
fMRI data acquisition 

fMRI data were acquired using a 3T GE Discovery MRI scanner 
with a T2* -weighted spiral BOLD sequence with pulse sequence 
parameters TR/TE/FA = 2 s/30 s/90°. The FOV was 22 cm in a 
64 x 64 matrix for 40 slices with 3 mm thickness. The total time 
for each scan (400 s) was matched to the functional task at 200 
volumes, after discarding 5 volumes at the beginning of each scan. 
Physiological signals were collected concurrently using a pulse 
oximeter and chest plethysmograph. For anatomical reference 
in the functional data analysis, a high-resolution Tl-weighted 
anatomical image was collected using spoiled-gradient-recalled 
acquisition (SPGR) in steady-state imaging with pulse sequence 
parameters TR/TE/FA = 12.2 ms/5.2 ms/15°. The FOV was 26 cm 
in a 256 x 256 matrix for 136 slices at thickness 1.2 mm. 

fMRI data processing 

After slice-timing and motion correction, and physiological noise 
removal with RETROICOR (Glover et al, 2000), fMRI data were 
co-registered to the participant's anatomical scan, normalized 
to the MNI template (Collins et al., 1994), smoothed, modeled 
and estimated with SPM8 (Wellcome Department of Cognitive 
Neurology, London, UK). Smoothing was performed with a 
Gaussian kernel of 8 x 8 x 8 mm at full width half maximum. 
The resulting image files contained BOLD activation time traces 
for each voxel of the brain and were used in the classification step. 

Both "rest minus task" and "task minus rest" contrasts were 
generated and used in second level random effects analyses per- 
formed across all runs for seven participants to inform fNIRS 
probe placement, and separately for each participant (across that 
individual participant's four runs) to guide voxel selection for 
fMRI classification. 

Across-network, within-network and co-activating correla- 
tions were defined, respectively, as the correlation coefficients 
between across-network (i.e., the DLPFC and the MFG), within- 
network (i.e., DLPFC and DLPFC2), or co-activating (i.e., the 
DLPFC and motor cortex) pairs of functional task signals. The 



Table 1 | The best c, g parameters found after optimization for each 
participant for fNIRS classification, and the time trace correlations (r) 
averaged across that participant's four runs. 



Participant 


c 


9 


r (within) 


r (across) 


1 


1 


0.1 


0.84 


-0.21 


2 


10 


0.001 


0.23 


-0.01 


3 


10 


0.01 


0.82 


-0.16 


4 


10 


0.1 


0.95 


-0.07 


5 


0.1 


0.001 


0.91 


-0.15 



fMRI traces were additionally smoothed across 10 s to reduce the 
impact of noise in the signal, before the correlation coefficient was 
calculated. These were averaged across runs by participant, then 
averaged across participants. 

fMRI classification 

All fMRI traces used as SVM inputs were processed BOLD 
responses, averaged across clusters centered on local maxima 
within the regions of interest (DLPFC, primary motor cor- 
tex, and MFG). Eighteen voxels were symmetrically selected 
around participant-specific centers to be included in the average. 
Participant-specific locations were selected using the second-level 
statistical maps generated using that individual participant's four 
runs. The voxels used were not contiguous, and the region of 
interest spanned 1 cm per side in MNI space. In this way, the 
SVM input features for fMRI were restricted to traces averaged 
across local tissue. This is analogous, for fairness of compar- 
ison between the modalities, to the volume of tissue interro- 
gated by one fNIRS probe, which is on the order of centimeters 
(Boas etal, 2004). 

Classification was performed with a fMRI trace selected from 
the contralateral motor area, which was paired with one from the 
DLPFC region (a co-activating pair, noted as DLPFC and Motor 
in Figure 4). The same DLPFC trace was paired with one from 
the MFG region (an across-network pair, noted as DLPFC and 
MFG in Figure 4), and classification was performed again. The 
same DLPFC trace was then paired with a second DLPFC trace 
from a region 2 cm superior and 2 cm medial to the first DLPFC 
trace in MNI space (a within-network pair, noted as DLPFC + 
DLPFC2 in Figure 4), and classification was performed a third 
time. To optimize accuracy, c was tested at 0.0001, 0.001, 0.01, 
0.1, 1, 5, 10,100,1000, and g was tested at 0.00001, 0.0001, 0.001, 
0.01,0.1,1, 5, 10, 100 (a total of 81 cases) in each classification step. 
Otherwise, all methods for classification were the same as those 
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FIGURE 4 | Classification accuracy averaged across four runs for each 
of seven participants. Two fMRI time traces were used as support vector 
machine input features in each case. Accuracy is the number of time points 
for which the prediction matched the truth label out of all 200 volumes. An 
accuracy of 50% represents prediction by chance. Error bars represent ± 
one standard deviation. Left-hand bars: co-activating pairs (blue). Center 
bars: within-network pairs (green). Right-hand bars: cross-region pairs 
(red). 
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Table 2 | The best c, g parameters found after optimization for each 
participant for fMRI classification, and the time trace correlations (r) 
averaged across that participant's four runs. 
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6 


5 


0.1 
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for the fNIRS traces described above. The SVM tuning parameters 
producing the best classification accuracies after optimization are 
given in Table 2. The accuracies are summarized in Figure 4. 

The motor region was selected for its robust and reliable 
response during task periods. We treat classification based on 
time traces from the motor region as a gold standard. That is, 
we do not expect fNIRS classification accuracies to exceed those 
attainable using motor cortex activations measured with fMRI. 

RESULTS 

BEHAVIORAL TASK RESULTS 

As described above, five participants completed the fNIRS study 
and seven participants completed the fMRI study. One partici- 
pant performed the task incorrectly for incongruent trials. Thus, 
this participant's data were excluded from the behavioral analy- 
ses. However, since this participant was engaged in the task and 
responding to stimuli, these trials were not excluded from the 
fNIRS and fMRI analyses. Including these data was appropriate 
because the fNIRS and fMRI analyses were aimed at distin- 
guishing between performing a task and resting, rather than 
distinguishing between incongruent and congruent trials. Mean 
accuracy was 98% (SD = 2.2%, N = 4) for the four fNIRS partic- 
ipants who performed the task correctly, and 97% (SD = 4.7%, 
N = 6) for the six fMRI participants who performed the task 
correctly. As expected, mean accuracy was relatively high. 

Also as expected, performance was worse in incongruent than 
in congruent trials (the analysis of the data included only the par- 
ticipants who performed correctly in most incongruent trials). In 
a random effects analysis, mean reaction time was significantly 
higher in the incongruent condition (fNIRS: M = 0.639 s, SD = 
0.035, N = 4; fMRI: M = 0.789 s, SD = 0.114 s, N = 6) than 
in the congruent condition (fNIRS: M = 0.552 s, SD = 0.033 s, 
N = 4; fMRI: M = 0.714 s, SD = 0.096 s, N = 6), [fNIRS: f (3) = 
26, p < 0.0005; fMRI: f (5) = 9.1, p < 0.0005]. Likewise, mean 
error rate was significantly higher in the incongruent condi- 
tion (fNIRS: M = 2.0%, SD = 1.8%, N = 4; fMRI: M = 4.9%, 
SD = 6.4%, N = 6) than in the congruent condition (fNIRS: 
M = 0.38%, SD = 0.81%, N = 4; fMRI: M = 2.1%, SD = 3.5%, 
N = 6), [fNIRS: f (3) = 2.9, p < 0.05; fMRI: f (5) = 2.2, p < 0.05). 
Errors of omission were rare and thus not analyzed. 

fNIRS RESULTS 

The four-run average classification accuracy for each of five par- 
ticipants is presented in Figure 2. As expected, we were able 



to distinguish between task engagement and rest. In particular, 
both averages differed significantly from chance at 50% [across: 
t(4) = 9.65, p < 0.0005; within: f (4) = 4.95, p < 0.005]. Also in 
line with predictions, there was a non-significant trend toward 
greater classification accuracy for across-network pairs (M = 
69.1%, SD = 4.4%) than for within-network pairs (M = 66.0%, 
SD = 7.2%), [f (4) = 1.48, p < 0.25]. The probe channels result- 
ing in oxygenated or deoxygenated [Hb] traces with the highest 
task model fit parameters were 1, 2 and 3 (as numbered in 
Figure 1, left; data not shown). The [Hb] species of the best- 
fitting traces were nearly evenly split: eight were oxygenated and 
seven were deoxygenated [Hb] traces. Thus, we found no uni- 
versally best Hb species for fitting the task model. This result is 
consistent with prior suggestions that a probe's sensitivity to one 
species or the other depends on whether it mostly samples the 
arterial or venous compartment (Strangman et al., 2003), which 
is likely to change for every probe application. 

Prior work indicates that activity in key regions of the 
"task-positive" network is positively correlated while activity in 
key regions of the "task-positive" and "task-negative" networks 
is negatively correlated (Fox et al, 2005; Kelly et al, 2008). 
Consistent with such findings, the group-averaged correlations 
for within-network pairs were significantly greater than zero [r = 
0.75, SD = 0.29; f (4) = 5.72, p < 0.005], while those for across- 
network pairs were significantly less than zero [r = —0.12, SD = 
0.08; f(4) = 3.41, p < 0.025]. Further, the correlation averages for 
within-network pairs were significantly higher than those of the 
across-network pairs [t^) = 5.56, p < 0.005; see Table 1 for a 
participant-specific list of correlation values] . These findings sug- 
gest that our fNIRS probes accurately measured activity in the 
"task-positive" and "task-negative" networks. 

Finally, since across-network correlation was determined on a 
per-participant basis, we also wished to verify that negatively cor- 
related across-network activity was observed at the group level. 
To this end, we averaged the across-network functional task sig- 
nals across all 20 runs. For the purposes of this group average, all 
deoxygenated traces were first inverted, consistent with a reduc- 
tion of deoxygenated [Hb] during activation. The group-averaged 
time traces are presented in Figure 3, which shows filtered, nor- 
malized and corrected [Hb] changes [see fNIRS methods; truth 
labels (green trace) show task at +1 and rest at —1]. As expected, 
the across-network correlation determined in this fixed effects 
analysis was significantly less than zero [r(i$) = -0.58, p < 0.01]. 
This finding illustrates that, even at the group level, DLPFC activ- 
ity (blue trace) increased during task performance while MFC 
activity (red trace) decreased. 

fMRI RESULTS 

The four-run average classification accuracy for each of the 
seven participants is presented in Figure 4. Replicating the fNIRS 
results, group averaged classification accuracy was significantly 
greater than chance at 50% for all three types of region pairs 
[across: f(6) = 3.56, p < 0.01; within: = 2.56, p < 0.025; co- 
activating: f(g) = 7.96, p < 0.0005], indicating we were able to 
distinguish task engagement from rest. Also as expected, group 
averaged classification accuracy was significantly higher for the 
co-activating (DLPFC and motor) pairs (M = 74.1%, SD = 
8.0%) than for the across-network pairs (M = 65.6%, SD = 
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11.6%), [f (6) = 4.83, p < 0.005], consistent with the high relia- 
bility of motor cortex activation detection in fMRI studies (Moller 
et al., 2005) and with the motor cortex activation associated with 
button press responses in the present task. Within-network pair 
accuracy (M = 63.0%, SD = 13.5%) tended to be lower than 
across-network pair accuracy, but not significantly [t($) = 1.57, 
p<0.1]. 

Also consistent with the fNIRS results, the group-averaged 
correlations for co-activating (r = 0.53, SD = 0.23) and within- 
network pairs (r = 0.54, SD = 0.26) were significantly greater 
than zero [within: f( 6 ) = 6.07, p < 0.0005; co-activating: f (6 ) = 
5.49, p < 0.005]. In contrast, those for the across-network pairs 
(r = 0.17, SD = 0.25) did not differ from zero [f (6) = 1.81, p < 
0.1]. As predicted, however, they were significantly lower than 
those for the within-network pairs [tm =2.29, p < 0.05; see 
Table 2 for a participant-specific list of correlation values]. 

The locations of statistically significant activations for the 
MSIT, after second-level analysis across four runs each from seven 
independent participants, are shown for the DLPFC [(52, 14, 32) 
in MNI space; Figure 5, left] and for the MFG [(22, 66, 0) in MNI 
space; Figure 5, right] . The f-statistic is mapped, with the thresh- 
old set at an uncorrected significance level of p < 0.001 for the 
work minus rest contrast shown on the left, and at p < 0.01 for 
the rest minus work contrast shown on the right. The MFG acti- 
vation was not present at the higher threshold but appeared at the 
lower threshold. Of importance, the expected "task-positive" and 
"task-negative" hemodynamic activations occurred in the same 
regions that were interrogated by the fNIRS probes. Notably, even 
with only seven participants, the DLPFC survived a family wise 
error correction at p < 0.05; the MFG, however, did not survive 
this correction. 

Finally, we note that although other regions of the "task- 
negative" network were identified in the fMRI analysis, they may 
be less useful for monitoring task engagement with fNIRS. First, 
lateral parietal regions were not consistently activated bilater- 
ally across participants. Thus, monitoring both sides with fNIRS 



would present greater difficulty due to the increased number of 
optical probes. Second, although precuneus and PCC regions of 
the task-negative network (Raichle et al., 2001; Greicius et al., 
2003) were reliably activated (Figure 5, right), they are too deep 
to be accessible via fNIRS probes, which can interrogate only the 
outer layers of the cortex (Boas et al., 2004). 

DISCUSSION 

In the present study, we investigated whether functional neu- 
roimaging methods (i.e., fNIRS and fMRI) can be employed to 
distinguish periods of task engagement from periods of rest. As 
described below, our findings support this view. They also provide 
valuable information about which brain activations may prove 
most useful for monitoring task engagement in the field. 

Our first set of findings came from fNIRS experiment. Here, 
we found that multivariate pattern classification techniques could 
distinguish between periods of task performance and periods of 
rest based on brain activity recorded from (a) different regions 
of the DLPFC in the task-positive network (a within-network 
pair) or (b) the DLPFC in the task-positive network and the MFG 
in the task-negative network (an across-network pair). Further, 
there was a trend toward higher classification accuracy for across- 
network pairs than for within-network pairs. Indeed, accuracy 
with across-network pairs approached 70%, even with the basic 
processing methods described here (with adaptive physiological 
filtering and additional probes, accuracy may further improve). 
This result fits with prior data suggesting that variability in task 
engagement is associated with variability in activity and/or func- 
tional connectivity involving both the "task-positive" and the 
"task-negative" networks (e.g., Weissman et al., 2006; Prado and 
Weissman, 2011). Most important, our fNIRS findings indicate 
that online recordings of brain activity via fNIRS may provide a 
valuable tool for detecting varying levels of task engagement in 
the real world. 

Our second set of findings came from an fMRI study. Of 
importance, these findings both verified and extended the results 



m 



FIGURE 5 | The locations of statistically significant activations for the 
MSIT, after second-level analysis across seven participants. The 

expected "task-positive" and "task-negative" hemodynamic activations 
occurred in the same regions that were interrogated by the fNIRS 
probes. Left: the f-statistic is mapped, with the height threshold set at 
an uncorrected significance level of p < 0.001 and the extent threshold 




set at 10 voxels. DLPFC is shown at [52, 14, 32] as marked by the 
crosshair in MNI space for the work minus rest contrast. Right: the 
f-statistic is mapped, with the height threshold set at an uncorrected 
significance level of p < 0.01 and the extent threshold set at 10 voxels. 
MFG is shown at [22, 66, 0] in MNI space for the rest minus work 
contrast. (Crosshair: x = 0, y = 66, z=14). 
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of the fNIRS study discussed earlier. First, we observed acti- 
vations and deactivations, respectively, in the DLPFC and the 
MFG, which verified that our functional neuroimaging paradigm 
engaged the task-positive and task-negative networks. Second, 
further analyses revealed that these activations occurred in the 
same DLPFC and MFG regions that were activated in the fNIRS 
experiment, wherein activity was measured with probes on the 
scalp. Third, the trend toward higher classification accuracy for 
across-network pairs than for within-network pairs observed with 
fNIRS was also observed with fMRI. Fourth, the fMRI findings 
indicated a possible upper bound on classification accuracy for 
distinguishing between task performance and rest: as expected, 
the highest classification accuracy was observed when the DLPFC 
trace and co-activating (for this task) motor cortex trace served as 
inputs to the SVM classifier. Together, these fMRI findings veri- 
fied that our fNIRS recordings reflected activity in the key regions 
under investigation (DLPFC and MFG). They also replicated the 
fNIRS results and extended them by suggesting an upper bound 
for classification accuracy based on relatively focal hemodynamic 
activity. 

NOVEL CONTRIBUTION OF THE PRESENT WORK 

The present findings make an important contribution to the field. 
Specifically, they show, for the first time, that it is possible to 
detect negative correlations between activity in key regions of the 
"task-positive" and "task-negative" networks with fNIRS. Further, 
they show that the detection of activity in the "task-negative" net- 
work is useful for distinguishing between high and low levels of 
task engagement. This capability might prove useful in future 
applications of fNIRS that are aimed at discriminating between 
optimal behavioral performance (where a negative correlation 
is expected) and internally-guided thought (where co-activation 
and, hence, a positive correlation is expected) (Christoff et al., 
2009; Smallwood et al., 2012). Thus, it could function to improve 
the predictive power of a fNIRS-based cognitive state monitoring 
system. Finally, although our findings make a novel contribu- 
tion to the field, it is important to note that they build on 
previous work showing that frontal oxygenation is sensitive to 
workload (Izzetoglu et al., 2004) and that fNIRS can reliably 
detect both resting state physiology and functionally-connected 
networks (White and Culver, 2008; Mehnert et al, 2009; Mesquita 
et al, 2010). 

LIMITATIONS 

While across-network pairs were associated with stronger neg- 
ative correlations and higher classification accuracy, relative to 
within-network measures (data not shown), negative correla- 
tions between DLPFC and MFG activity were not observed in 
every participant. This lack of consistency may stem from a 
variety of sources, including non-optimal fNIRS probe localiza- 
tion, variation in participant compliance or strategy, interference 
from physiological or motion artifact, variable fMRI voxel selec- 
tion, and co-activation of key regions in the "task-positive" and 
"task-negative" networks during mind wandering or internally- 
guided thought (Christoff et al, 2009; Smallwood et al., 2012). 
Future studies should be conducted to distinguish among these 
possibilities and to determine which methodologies provide more 



consistent measures of negatively correlated activity in the task- 
positive and task-negative networks. 

Also regarding the consistency of our measures, classification 
accuracy varied considerably across runs (see the error bars in 
Figures 2, 4). Future studies might therefore be conducted to 
investigate the source(s) of this variability as well as the impact of 
other sources of variability (e.g., across-visit, across-participant 
and across-task) on classification accuracy. Such studies might 
also investigate the impact of using the known task model to clean 
the measured traces when producing functional task signals for 
use in classification (see fNIRS data processing), which may have 
biased the classifier toward higher accuracy in the present study. 

Another study limitation stems from the fact that some task- 
evoked systemic signals are measureable on the scalp surface, and 
that at least one such signal — skin blood volume — depends on 
cognitive state (Kirilina et al., 2012). Since fNIRS is sensitive to 
hemodynamics in superficial tissue at all source-detector separa- 
tion distances, it is possible that systemic signals in the superficial 
tissue may have driven the negative across-network correlations 
that we observed. To investigate this possibility, we quantified cor- 
relations for the traces taken from the across-network shallow 
source-detector pairs (see fNIRS data acquisition). Of impor- 
tance, no association was observed between the correlation values 
for the superficial traces and the correlation values for the deep 
traces, whether corrected or not (data not shown). Thus, the neg- 
ative correlations that we measured between the across-network 
deep traces likely reflected [Hb] related changes in brain tis- 
sue rather than superficial physiological signals. However, if skin 
blood changes provide additional information about the task 
engagement, then it could be useful in future studies to include 
such changes directly as classifier inputs. 

Finally, we note that motor activation is not always a reliable 
component of task engagement as some tasks require sustained 
attention over long periods in the absence of overt responses (e.g., 
instrument cross-checking and visual display searches). Thus, 
future studies may wish to focus on our fNIRS finding that classi- 
fication accuracy was slightly higher for across-network pairs (i.e., 
pairs in which one region came from the task-positive network 
while the other came from the task-negative network) than for 
within-network pairs (i.e., pairs in which both regions came from 
the task-positive network). As we mentioned earlier, this finding 
fits with the view that task engagement is determined by interac- 
tions between these networks (Fox et al, 2005; Weissman et al., 
2006; Kelly et al, 2008; Prado and Weissman, 2011). 

FUTURE WORK 

Future work could address whether adaptive filtering of fNIRS 
traces over smaller time windows improves the activation-based 
classification measures reported here. This may include investi- 
gating adaptive physiological noise removal driven by the correla- 
tion between the deep and shallow traces (Harrivel et al., 2012) or 
motion artifact reduction based on the frequency domain phase 
signal (Harrivel and Hearn, 2012). An increase in probe density 
could also be used to improve localization within the regions of 
interest. Finally, measures of network correlation could be quan- 
tified over shorter time scales to determine whether transient 
internally- guided thought can be distinguished from periods of 
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"zoning out." Ongoing simultaneous fNIRS/fMRI studies are 
further examining these and other possible methods for improv- 
ing our ability to discriminate between varying levels of task 
engagement. 

CONCLUSION 

In the present study, we used a combination of fNIRS and fMRI 
results to show that online recordings of brain activity from the 
task-positive and task-negative networks can be used to detect 
moment-to-moment changes in task engagement. We hope that 
future studies combining fNIRS recordings with multivariate clas- 
sification methodologies will build upon the present work to 
produce robust systems for detecting changes in task engagement 
in real-world settings. 

AUTHORS CONTRIBUTIONS 

Angela R. Harrivel is the PI on the human subject study. 
She recruited and consented participants, designed fNIRS head 
probes, wrote all scripts for the fNIRS data processing, col- 
lected and processed all data, performed analyses and wrote the 
manuscript. Daniel H. Weissman wrote the code for the MSIT 
presentation, provided guidance regarding psychological aspects 
of the study and fMRI analyses, and edited the manuscript. 
Douglas C. Noll provided funding support, guidance regarding 
all fMRI aspects of the project and all data processing methods, 
and reviewed the manuscript. Scott J. Peltier is Co-I on the study. 
He contributed to analyses, edited the manuscript, and provided 
guidance regarding data collection and processing techniques, 
pattern classification methods, and resting state analyses. 

ACKNOWLEDGMENTS 

This work was supported by the University of Michigan fMRI 
Laboratory and NASA's Aviation Safety Program. Colleagues at 
the NASA Glenn and Langley Research Centers are appreci- 
ated, especially leffrey Mackey, Daniel Gotti and Padetha Tin for 
head probe design and assembly, and Tristan Hearn and Alan 
Pope for helpful review. We are grateful for the assistance of the 
fMRI laboratory at the University of Michigan for the collec- 
tion and pre-processing of the fMRI data, and Ted Huppert of 
the University of Pittsburgh for review and invaluable guidance 
regarding general fNIRS techniques and analyses. 

REFERENCES 

Boas, D. A., Dale, A. M., and Franceschini, M. A. (2004). Diffuse optical imaging 

of brain activation: approaches to optimizing image sensitivity, resolution, and 

accuracy. Neuroimage 23, S275-S288. doi: 10.1016/j.neuroimage.2004.07.011 
Bunce, S., Izzetoglu, M., Izzetoglu, K., Onaral, B., and Pourrezaei, K. (2006). 

Functional near-infrared spectroscopy — an emerging neuroimaging modality. 

IEEE Eng. Med.Biol Mag. 25, 54-62. doi: 10.1109/MEMB.2006.1657788 
Bush, G., and Shin, L. (2006). The multi-source interference task: an fMRI task that 

reliably activates the cingulo-frontal-parietal cognitive/attention network. Nat. 

Protoc. 1, 308-313. doi: 10.1038/nprot.2006.48 
Chang, C. C, and Lin, C. J. (2011). LIBSVM: a library for support vector machines. 

ACM Trans. Intel!. Syst. Technol 2, 1-27. doi: 10.1145/1961189.1961199 
Chee, M. W. L., Tan, J. C, Zheng, H., Parimal, S., Weissman, D. H., 

Zagorodnov, V., et al. (2008). Lapsing during sleep deprivation is associated 

with distributed changes in brain activation. /. Neurosci. 28, 5519-5528. doi: 

10.1523/JNEUROSCI.0733-08.2008 
Christoff, K., Gordonb, A. M., Smallwood, J., Smith, R., and Schooler, J. W. 

(2009). Experience sampling during fMRI reveals default network and exec- 



utive system contributions to mind wandering. PNAS. 106, 8719-8724. doi: 
10.1073/pnas.0900234106 
Collins, D. L., Neelin, P., Peters, T. M., and Evans, A. C. (1994). Automatic 3D 
intersubject registration of MR volumetric data in standardized Talairach space. 
/. Comput. Assist. Tomogr. 18, 192-205. doi: 10.1097/00004728-199403000- 
00005 

Cowings, P. S., Toscano, W. B., Taylor, B., DeRoshia, C. W., Kornilova, L. 
Kozlovskya, I., et al. (2003). "Psychophysiology of spaceflight," 14th IAA 
Humans In Space Symposium, Living in Space: Scientific, Medical and Cultural 
Implications (Banff, Alberta). 

De Havas, J. A., Parimal, S., Soon, C. S., and Chee, M. W. L. (2012). 
Sleep deprivation reduces default mode network connectivity and anti- 
correlation during rest and task performance. Neuroimage 59, 1745-1751. doi: 
10.1016/j.neuroimage.201 1.08.026 

Delpy, D. T, Cope, M., van der Zee, P., Arridge, S., Wray, S., and Wyatt, J. (1998). 
Estimation of optical pathlength through tissue from direct time of flight mea- 
surement. Phys. Med. Biol 33, 1433-1442. doi: 10.1088/0031-9155/33/12/008 

Dismukes, R. K., Berman, B. A., and Loukopoulos, L. D. (2007). The Limits of 
Expertise: Rethinking Pilot Error and the Causes of Airline Accidents (Ashgate 
Studies in Human Factors for Flight Operations). Burlington, VT: Ashgate 
Publishing Company. 

Dosenbach, N. U. E, Visscher, K. M., Palmer, E. D., Miezin, F. M., Wenger, K. K., 
Kang, H. C, et al. (2006). A core system for the implementation of task sets. 
Neuron 50, 799-812. doi: 10.1016/j.neuron.2006.04.031 

Drummond, S.P. A., Bischoff-Grethe, A., Dinges, D. E, Ayalon, L., Mednick, S. C, 
and Meloy, M. J. (2005). The neural basis of the psychomotor vigilance task. 
Sleep. 28, 1059-1068. 

Emir, U. E., Ozturk, C. and Akin, A. (2008). Multimodal investigation of fMRI 
and fNIRS derived breath hold BOLD signals with an expanded balloon model. 
Physiol Meas. 29, 49-63. doi: 10.1088/0967-3334/29/1/004 

Fox, M. D., Snyder, A. Z., Vincent, J. L., Corbetta, M., Van Essen, D. C, 
and Raichle, M. E. (2005). The human brain is intrinsically organized into 
dynamic, anticorrelated functional networks. PNAS 102, 9673-9678. doi: 
10.1073/pnas.0504136102 

Gagnon, L., Perdue, K., Greve, D. N., Goldenholz, D., Kaskhedikar, G., and Boas, 
D. (2011). Improved recovery of the hemodynamic response in diffuse optical 
imaging using short optode separations and state-space modeling. Neuroimage 
56, 1362-1371. doi: 10.1016/j.neuroimage.201 1.03.001 

Gibson, A. P., Hebden, J. C, and Arridge, S. R. (2005). Recent advances in diffuse 
optical imaging. Phys. Med. Biol 50, R1-R43. doi: 10.1088/0031-9155/50/4/R01 

Glover, G. H., Li, T. Q., and Ress, D. (2000). Image-based method for retrospec- 
tive correction of physiological motion effects in fMRI: RETROICOR. Magn. 
Reson. Med. 44, 162-167. doi: 10.1002/1522-2594(200007)44:1<162::AID- 
MRM23>3.3.CO;2-5 

Gratton, E., Toronov, V, Wolf, U, Wolf, M., and Webb, A. (2005). Measurement 
of brain activity by near-infrared light. /. Biomed. Opt. 10:011008-1-13. doi: 
10.1117/1.1854673 

Greicius, M. D., Krasnow, B., Reiss, A. L., and Menon, V. (2003). Functional con- 
nectivity in the resting brain: a network analysis of the default mode hypothesis. 
Proc. Natl Acad. Sci. U.S.A. 100:1, 253-258. doi: 10.1073/pnas.0135058100 

Harrivel, A., and Hearn, T. (2012). "Functional near infra-red spectroscopy: watch- 
ing the brain in flight," in 4th International Conference on Applied Human Factors 
and Ergonomics (San Francisco, CA). 

Harrivel, A., Hernandez-Garcia, L., Peltier, S., and Noll, D. (2012). "Artifact 
removal for assessment of cross-network anticorrelation and attentional clas- 
sification with fNIRS," in 2nd Biennial fNIRS Meeting (London). 

Huppert, T. J., Diamond, S. G., Franceschini, M. A., and Boas, D. A. (2009). 
HomER: a review of time-series analysis methods for near-infrared spec- 
troscopy of the brain. Appl. Opt. 48, D280-D298. doi: 10.1364/AO.48.00D280 

Huppert, T. J., Hoge, R. D., Dale, A. M., Franceschini, M. A., and Boas, D. A. (2006). 
Quantitative spatial comparison of diffuse optical imaging with blood oxygen 
level-dependent and arterial spin labeling-based functional magnetic resonance 
imaging./. Biomed. Opt. 11, D280-D298. doi: 10.1117/1.2400910 

Izzetoglu, K., Bunce, S., Onaral, B., Pourrezaei, K., and Chance, B. (2004). 
Functional optical brain imaging using near-infrared during cognitive tasks. Int. 
}. Hum.Comput. Interact. 17, 211-231. doi: 10.1207/sl5327590ijhcl702_6 

Jimenez, L., and Mendez, A. (2013). It is not what you expect: dissociating conflict 
adaptation from expectancies in a Stroop task. /. Exp. Psychol. Hum. Percept. 
Perform.39, 271-284. doi: 10.1037/a0027734 



Frontiers in Human Neuroscience 



www.frontiersin.org 



December 2013 | Volume 7 | Article 861 | 9 



Harrivel et al. 



Monitoring attentional state with fNIRS 



Kelly, A. M. C, Uddin, L. Q., Biswal, B. B., Castellanos, F. X., and Milham, M. P. 
(2008). Competition between functional brain networks mediates behavioral 
variability. Neuroimage 39, 527-537. doi: 10.1016/j.neuroimage.2007.08.008 

Kirilina, E., Jelzowb, A., Heinea, A., Niessinga, M., Wabnitzb, H., Briihlb, 
R., et al. (2012). The physiological origin of task-evoked systemic arte- 
facts in functional near infrared spectroscopy. Neuroimage 61, 70-81. doi: 
10.1016/j.neuroimage.2012.02.074 

Kleiner, M., Brainard, D., and Pelli, D. (2007). What's new in psychtoolbox-3? 
Perception 36. (ECVP abstract supplement), 14. doi: 10.1068/v070821 

Kleinschmidt, A, Obrig, H., Requardt, M., Merboldt, K, Dirnagl, U., Villringer, 
A., et al. (1996). Simultaneous recording of cerebral blood oxygenation 
changes during human brain activation by magnetic resonance imaging and 
near-infrared spectroscopy. /. Cereb. Blood Flow Metab. 16, 817-826. doi: 
1 0. 1097/00004647- 1 99609000-00006 

Li, L., Du, R, Li, T., Luo, Q., and Gong, H. (2007). "Design and evaluation of 
a simultaneous fNIRS/ERP instrument," in Conference on Optical Tomography 
and Spectroscopy of Tissue VII. Proceedings of SPIE. 6434 (643429) (San Jose, 
CA). 

MacDonald, A. W., Cohen, J. D., Stenger, V. A, and Carter, C. S. (2000). 
Dissociating the role of the dorsolateral prefrontal and anterior cingu- 
late cortex in cognitive control. Science 288, 1835-1838. doi: 10. 1 126/sci- 
ence.288.5472.1835 

MATLAB. (2012). Statistics and Signal Processing Toolboxes, R2012b. Natick, MA: 
The Math Works, Inc. 

Mayr, U., Awh, E., and Laurey, P. (2003). Conflict adaptation effects in the absence 
of executive control. Nat. Neurosci. 6, 450-452. doi: 10.1038/nnl051 

McKiernan, K. A. Kaufman, J. N., Kucera-Thompson, J., and Binder, J. R. 
(2003). A parametric manipulation of factors affecting task-induced deac- 
tivation in functional neuroimaging. /. Cogn. Neurosci. 15, 394-408. doi: 
10.1162/08989290332159317 

Mehnert, J., Margulies, D. S., Schmitz, C, Steinbrink, J., Obrig, H., and Villringer, 
A. (2009). "Resting state networks revealed with whole-head near-infrared 
spectroscopy," in I5th Annual Meeting of the Organization for Human Brain 
Mapping (San Francisco, CA). Neuroimage 47, S163. doi: 10.1016/S1053-8119 
(09)71728-8 

Mesquita, R. C, Franceschini, M. A., and Boas, D. A. (2010). Resting state func- 
tional connectivity of the whole head with near-infrared spectroscopy. Biomed. 
Opt. Express, 1, 324-336. doi: 10.1364/BOE.1.000324 

Moller, M., Freund, M., Greiner, C, Schwindt, W., Gaus, C, and Heindel, W. 
(2005). Real time fMRI: a tool for the routine presurgical localisation of the 
motor cortex. Eur. Radiol. 15, 292-295. doi: 10.1007/s00330-004-2513-z 

Moosmann, M., Ritter, P., Krastel, I., Brink, A., Thees, S., Blankenburg, E, et al. 
(2003). Correlates of alpha rhythm in functional magnetic resonance imaging 
and near infrared spectroscopy. Neuroimage 20, 145-158. doi: 10.1016/S1053- 
8119(03)00344-6 

Obrig, EL, and Villringer, A. (2003). Beyond the visible — imaging the human 
brain with light. /. Cereb. Blood Flow Metab. 23, 1-18. doi: 10.1097/00004647- 
200301000-00001 

Pope, A. T, Bogart, E. H., and Bartolome, E. S. (1995). Biocybernetic system 

evaluates indices of operator engagement in automated task. Biol. Psychol. 40, 

187-195. doi: 10.1016/0301-0511(95)05116-3 
Prado, J., and Weissman, D. H. (2011). Heightened interactions between a key 

default-mode region and a key task-positive region are linked to suboptimal 

current performance but to enhanced future performance. Neuroimage 56, 

2276-2282. doi: 10.1016/j.neuroimage.2011.03.048 
Raichle, M. E., MacLeod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and 

Shulman, G. L. (2001). A default mode of brain function. PNAS 98, 676-682. 

doi: 10.1073/pnas.98.2.676 
Raley, C, Stripling, R., Kruse, A., Schmorrow, D, and Patrey, J. (2004). 

Augmented cognition overview: improving information intake under 



stress. Proc. Hum. Factors Ergon. Soc. Annu. Meet. 48, 1150-1154. doi: 
10.1177/154193120404801001 

Schnell, T., Kwon, Y., Merchant, S., and Etherington, T. (2004). Improved 
flight technical performance in flight decks equipped with synthetic 
vision information system displays. Int. J. Aviat. Psychol. 14, 79-102. doi: 
10.1207/sl5327108ijapl401_5 

Schroeter, M. L., Kupka, T, Mildner, T., Uludag, K, and Yves von Cramon, 
D. (2006). Investigating the post-stimulus undershoot of the BOLD signal — 
a simultaneous fMRI and fNIRS study. Neuroimage 30, 349-358. doi: 
10.1016/j.neuroimage.2005.09.048 

Smallwood, J., Brown, K, Baird, B., and Schooler, J. W. (2012). Cooperation 
between the default mode network and the frontal-parietal network in the 
production of an internal train of thought. Brain Res. 1428, 60-70. doi: 
10.1016/j.brainres.2011.03.072 

Steinbrink, J., Villringer, A., Kempf, F., Haux, D., Boden, S., and Obrig, 
H. (2005). "Illuminating the BOLD signal: combined fMRI-fNIRS stud- 
ies," in International School on Magnetic Resonance and Brain Function 
(Erice). 

Stins, J. E, van Leeuwen, W. M. A., and de Geus, E. J. C. (2005). The multi-source 

interference task: the effect of randomization. /. Clin. Exp. Neuropsychol. 27, 

711-717. doi: 10.1080/13803390490918516 
Strangman, G., Boas, D. A., and Sutton, J. P. (2002a). Non-invasive neuroimaging 

using near-infrared light. Biol. Psychiatry. 52, 679-693. doi: 10.1016/S0006- 

3223(02)01550-0 

Strangman, G, Culver, J. P., Thompson, J. EL, and Boas, D. A. (2002b). A 
quantitative comparison of simultaneous BOLD fMRI and NIRS record- 
ings during functional brain activation. Neuroimage 17, 719-731. doi: 
10.1006/nimg.2002.1227 

Strangman, G, Franceschini, M. A., and Boas, D. A. (2003). Factors affect- 
ing the accuracy of near-infrared spectroscopy concentration calculations 
for focal changes in oxygenation parameters. Neuroimage 18, 865-879. doi: 
10.1016/S1053-81 19(03)00021-1 

Strayer, D. L., Watson, J. M., and Drews, F. A. (2011). "Cognitive distraction while 
multitasking in the automobile," in The Psychology of Learning and Motivation, 
Vol. 54. ed B. Ross (Burlington, VT: Academic Press), 29-58. 

Weissman, D. FL, Roberts, K. C, Visscher, K. M. and Woldorff, M. G. (2006). The 
neural bases of momentary lapses in attention. Nat. Neurosci. 9, 971-978. doi: 
10.1038/nnl727 

White, B., and Culver, J. (2008). "Functional connectivity in adult humans revealed 
with diffuse optical tomography of oxy-, deoxy-, and total hemoglobin," in 14th 
Annual Meeting of the Organization for Human Brain Mapping (Melbourne, 
VIC). 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 30 September 2013; paper pending published: 22 October 2013; accepted: 26 
November 2013; published online: 13 December 2013. 

Citation: Harrivel AR, Weissman DH, Noll DC and Peltier SJ (2013) Monitoring 
attentional state with fNIRS. Front. Hum. Neurosci. 7:861. doi: 10.3389/fnhum. 
2013.00861 

This article was submitted to the journal Frontiers in Human Neuroscience. 
Copyright © 2013 Harrivel, Weissman, Noll and Peltier. This is an open-access article 
distributed under the terms of the Creative Commons Attribution License (CC BY). 
The use, distribution or reproduction in other forums is permitted, provided the 
original author(s) or licensor are credited and that the original publication in this 
journal is cited, in accordance with accepted academic practice. No use, distribution or 
reproduction is permitted which does not comply with these terms. 



Frontiers in Human Neuroscience 



www.frontiersin.org 



December 2013 | Volume 7 | Article 861 | 10 



